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Abstract 

The widely used nonperturbative wave functions and distribution functions 
of QCD are determined as matrix elements of light-ray operators. These op- 
erators appear as large momentum limit of nonlocal hadron operators or as 
summed up local operators in light-cone expansions. Nonforward one-particle 
matrix elements of such operators lead to new distribution amplitudes describ- 
ing both hadrons simultaneously. These distribution functions depend besides 
other variables on two scaling variables. They are applied for the description 
of exclusive virtual Compton scattering in the Bjorken region near forward 
direction and the two meson production process. The evolution equations 
for these distribution amplitudes are derived on the basis of the renormaliza- 
tion group equation of the considered operators. This includes that also the 
evolution kernels follow from the anomalous dimensions of these operators. 
Relations between different evolution kernels (especially the Altarelli-Parisi 
and the Brodsky-Lepage) kernels are derived and explicitly checked for the 
existing two-loop calculations of QCD. Technical basis of these results are 
support and analytically properties of the anomalous dimensions of light-ray 
operators obtained with the help of the a-representation of Green's functions. 
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I. INTRODUCTION 



It is generally accepted that quantum chromodynamics (QCD) allows the theoretical 
description of scattering processes at large momentum transfer. Thereby, it is possible to 
calculate the so-called hard scattering subprocess perturbatively, and for the remaining parts 
of the process one uses "parton distribution functions" |IJ or "wave functions" 0]. These 
functions are not calculable perturbatively, but have to be determined phenomenologically. 
Nevertheless additional information on them can be obtained because they satisfy evolution 
equations. 

In this paper we take into account that the usually applied nonperturbative distribution 
and wave functions are directly related to special matrix elements of (nonlocal) light-ray 
operators [0-0]. The parton distribution amplitude (corresponding to the forward matrix 
element) is necessary for the description of inclusive lepton-hadron scattering and the meson 
(hadron) wave functions (corresponding to the matrix element between the one-hadron and 
the vacuum state) is needed for the treatment of exclusive large transverse momentum pro- 
cesses. Here we discuss the question what role plays the general one-particle matrix element 
between hadrons of different momenta. To these matrix elements belong new distribution 
amplitudes describing both hadrons simultaneously. Such nonperturbative distribution am- 
plitudes are necessary if there is no large momentum transfer between the hadrons (otherwise 
we could use a a hard scattering part and a separate description of each hadron by a non- 
perturbative function). As an example we show that these functions can be applied to the 
prescription of the virtual Compton scattering amplitude in the Bjorken region near forward 
scattering or to two meson production processes. Thereby it is extremely interesting that 
we obtain a very simple representation of the Compton scattering amplitude in leading or- 
der. In fact, the virtual Compton amplitude can be represented as a convolution of a hard 
scattering part with the distribution amplitude, so that the analogy to the simple parton 
picture remains. 

A second topic of this paper is the investigation of the evolution equations and evolution 
kernels for these nonperturbative distribution amplitudes. Here we take into account that 
all these amplitudes are matrix elements of light-ray operators, satisfying renormalization 
group equations. Therefore, the evolution equations are basically renormalization group 
equations. Of course, by forming matrix elements the renormalization group equations for 
light-ray operators turn over to evolution equations for these matrix elements. Using such 
a general procedure, we are able (in principle) to derive evolution equations for all matrix 
elements of such operators. The Brodsky-Lepage equation [Q] for the meson wave function 
and the Altarelli-Parisi equation [[3J for the quark distribution function of deep inelastic 
scattering appear as special cases. It is important that the evolution kernels are functionals 
of the anomalous dimension of the considered operator. A trivial but interesting point is 
that thereby the range of the variables may change according to the support properties of 
the distribution amplitude. A very general kernel which is in fact the complete anomalous 
dimension of the underlying light-ray operator is the extended Brodsky-Lepage (BL) kernel. 
It represents a partly diagonalized kernel generating all one-parametric evolution equations. 
In the one-loop approximation, a complete diagonalized form of such a kernel following from 
nonlocal (conformal) operator product expansions has been studied in Ref. 0. 

Besides of these general considerations we discuss also some interesting technical prob- 
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lems concerning the explicit determination of these kernels. In QCD, the BL-kernel is calcu- 
lated as an evolution kernel in a restricted range of the variables. The interesting question 
is: Exists there a possibility to determine the extended BL-kernel knowing this restricted 
kernel only? We prove that the partly Fourier transformed BL-kernel has holomorphic prop- 
erties allowing such an extension. In the following these results are applied to the one- and 
two-loop calculations of the BL-kernel in QCD. Without additional diagrammatic calcula- 
tions we are able to determine the extended BL-kernel |7j] of QCD. Finally, we perform the 
limiting procedure leading from the extended BL-kernel to the Altarelli-Parisi (AP) kernel, 
and compare the resulting expression with the directly calculated results p-pH|]. In this way 
we obtain a consistency check for the very complicated two-loop calculations 0|1OJ. Taking 



into consideration the recent check [11] based on conformal Wardidentities and consistence 



equations, we conclude that the BL-kernel has been correctly computed. 

All our investigations concern the flavour nonsinglet operators. Investigations of singlet 
operators are in principle possible, however technically more complicated. 

It remains to give some remarks concerning the light-ray operators and its anomalous 



dimensions. Light-ray operators were first introduced by O. I. Zavialov [JT2J] in connection 
with his proof of the light-cone operator product expansion. In opposite to the standard 
formulation []T3| he did not use the expansion in terms of local operators |T3[ but introduced 
light-ray operators which are in fact summed up local operators. We also use this technique 
for the treatment of the virtual Compton scattering amplitude. We hope to convince the 
reader that this method is very effective. An independent approach based on the same 
ideas was developed and applied from the Leningrad group Another source for the 

introduction of light-ray operators are the speculations about hadron operators in QCD 



15| . Hadrons are composed of quarks and gluons, therefore hadron operators have to be bi- 
or trilocal operators containing the basic fields of QCD. The postulate of gauge invariance 
leads to string operators finally. If we want to apply these operators in hard scattering 
processes then for large internal momenta these operators turn over to the light-ray operators 
considered here. It is now clear that for example bilocal operators where the distance between 
both positions of the involved field operators is light-like possess additional divergences. 
Therefore they need a special renormalization procedure which leads to more complicated 
anomalous dimensions as usually discussed. 

To avoid complicated notations in the main text we generally write down renormalized 
expressions. Only in the Appendix we also use unrenormalized expressions. An important 
technical task is the investigation of the anomalous dimensions of these operators, especially 
the correct determination of its domain of definition. Using the a-representation of Feynman 
graphs, we first study this problem for the scalar theory and afterwards for QCD in axial 
gauge. For the last case we modify the a-representation to include the more complicated 
gluon propagator of QCD in light-cone gauge. 

The paper is organized as follows: In the second section we introduce the generalized 
distribution amplitudes as matrix elements of light-ray operators. In the next section we 
apply these distribution amplitudes for the prescription of the exclusive virtual Compton 
scattering and the two meson production by photons. In the fourth section we derive the 
evolution equations for these amplitudes. Problems concerning explicit expressions (two-loop 
calculations) and relation between different evolution kernels are discussed in the last section. 
Important technical details concerning the domain of definition of the anomalous dimension, 
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the perturbative calculation of hard scattering parts and other problems are shifted to 
Appendices. The Appendices are needed for the proofs only, the reader uninterested in 
technical details could skip them. 

II. NONPERTURBATIVE INPUTS: DISTRIBUTION AMPLITUDES 

Up to now the application of nonperturbatively determined "parton distribution func- 
tions" and "wave functions" is an unavoidable part of perturbative QCD calculations. We 
start our considerations with the introduction of more general distribution amplitudes cor- 
responding to nonforward scattering processes. 

In the QCD picture hadrons are built up from quarks and gluons. So possible hadron 
operators are nonlocal operators containing quark and gluon fields. Typical meson operators 
are composed of quark and antiquark operators 

O a (x 1 ,x 2 ) = ■4{x l )TX a U{x u x 2 )^{x 2 ): . (2.1) 

The path ordered phase factor, 

U(x 1 ,x 2 ) = Vexp |-z 5 p dx^ix)^ (2.2) 

ensures the gauge invariance of the considered operator. = A^t c denotes the gluon field, 
t c are the generators of the colour group, if) is the quark field, T symbolizes the necessary 
spin structure and A a is a generator of the flavour group corresponding to the considered 
meson. At large internal quark momenta the operator ( |2.1| ) turns over to an asymptotic 
operator defined on a straightline light-like path 

O a (n + , n) = : ^(Kin)(n7)A a f/(/tifi, K, 2 n)if)(K, 2 n) :, «± = (k 2 ± «i)/2, (2.3) 

where 

U(nih, K 2 n) = Pexp j— ig J drnA(rn)j. 

The vector n, with n 2 = 0, defines the light-ray pointing into the direction of the large 
momentum flow of this process. 

For the description of mesons involved in hard QCD processes one needs an "asymptotic" 
wave function as nonperturbative input. This meson wave function $ a (a;, Q 2 ) depending on 
the distribution parameter x (0 < x < 1) and the momentum transfer Q 2 can be defined as 
the expectation value of a nonlocal (light-ray) operator lying on the light-cone 0-|4|]: 

/n k \7i i 1 
' l e ^-(nP)tJ_ < Ql au_. ^ip > I a=g2 (2 .4) 
2n nP 

where 

O a ( K _; n) = O a (n + = 0, h) = :ip(-K_n)(n-f)\ a U(-K-fi, K^fi)ip(n„fi) : . (2.5) 

Here \P > denotes the one-particle state of a scalar meson of momentum P and O a (K_; n) is 
the light-ray operator with the same flavour content. As renormalization point we choose the 
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typical large momentum Q of the basic process to which the wave function contributes. The 
factor l/(nP) was introduced to compensate the n-dependence of the factor fry in the 
operator O a (K^;n). Contrary to the definition in Ref. our expression ( |2.4| ) is gauge 
invariant and Lorentz-covariant. In an infinite momentum frame with P = (\P Z \, 0, 0, P z ) 
(neglecting the masses of the particles) and h = (1,0,0,-1) the distribution parameter 
x can be interpreted as usually: xP is the fraction of the meson momentum P which is 
carried by the quark. Then the mesonwave function is the probability amplitude for finding 
a quark-antiquark pair in the meson in dependence of the momentum fraction parameter x 
and the momentum transfer Q 2 . 

As next, we define the quark distribution functions q^z, Q 2 ) [1[| with the distribution pa- 
rameter z. These functions are necessary for the description of deep inelastic scattering. At 
first such functions have been introduced phenomenologically, later on as Mellin transforms 
of moments of forward matrix elements of local operators [|IJ. If we take into account that 
the operator 

Oj(/« + ,K_) = \^ i (Kih)(n^)U(Kin,K2n)^ i (K 2 n): (2.6) 
is the "generating operator" for all local operators, 

{{nDy^m){ni){nDy^m = ——0^+, K _)| Kl =. 2 =o, (2.7) 

then a definition of the quark distribution function as Fourier transform of the one-particle 
matrix element of the light-ray operator 0], 

Oi(K_;n) = \^i(—K_h)(n^)U(—K_h 1 K_n)%p i (K_h)\ 1 (2-8) 

seems to be understandable. Here ipi denotes the quark field with flavour component i. 
However, group theoretically this operator has no definite flavour content which leads to 
mixing problems with the flavour singlet gluon operator. To avoid this problem we consider 
here the flavour nonsinglet distribution function only: 

q a {z,Q 2 ) = f dK -\ nP \ e ^-(nP)zJ_ < P \o a { K -;h)\P > |m 2 =Q 2 - (2.9) 

Physically this distribution function describes the probability of finding internal quarks 
with the momentum fraction zP, where P is the external hadron momentum. If the matrix 
A a is choosen diagonal, then for positive distribution parameter z this function represents a 
linear combination of the quark distribution functions, and the extension for z < can be 
interpreted as a linear combination of the antiquark distribution functions. 

Both functions in ( |2.4[ ) and ( |2.9| ) are expectation values of the same operator used, 
however, for different physical processes and having different interpretations. 

If we look at the definition (23!) of the quark distribution function given as Fourier 



transform of a forward matrix element of a light-ray operator then the question arises: 
What is the physical meaning of nonforward matrix elements? As an example we study in 
the next section two-photon processes in the Bjorken region. The idea is that the product 
of the two electromagnetic currents near light-like distances can be expressed with help of 
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coefficient functions and the light-ray operators O a (K_;x). (This generalizes the standard 
local operator product expansion). 

If we consider the virtual Compton scattering amplitude in nonforward direction, than we 
have to form nonforward matrix elements of these light-ray operators (see the next section) . 
This leads us directly to the following distribution amplitude, 

<f (t,r,/i 2 ) = J dK -\^ P +\ e ^-(nP + )t^_ < p 2 | a (K _. ~ )|pi > |- p _ =T - p+) (2.10) 

with P± = P 2 ± P\. This function depends on the distribution parameter t and the quo- 
tient r = nP_/nP + of the projection of momenta P± onto a light-like direction n, on the 
renormalization point \i and the scalar products of the external momenta PjPj. For physical 
states | Pi > and |P 2 > the additional variable r is restricted by 



T 







p o 




nP+ 




n 


-pi 



„ np , 

< 1, P" = — . (2.11) 

Inl 



With the help of the a-representation it is possible to investigate the support properties 
of the function q a (t,T,fi 2 ) with respect to the variable t. Extending the proof given in Ref. 
0j for t = 1 to arbitrary values of r, it turns out that 

q a (t,T,Q 2 ) = for \t\ > 1. (2.12) 

The distribution amplitude q a (t,r,Q 2 ) is especially suited for the description of nonfor- 
ward processes near to the forward case. Namely, if there is no large momentum transfer 
between the two hadrons, then it is impossible to describe the hadrons by individual wave 
functions and a hard scattering part, i.e. the parton picture with two separate wave func- 
tions cannot be applied. Here one has to apply one function describing both hadrons. An 
example of such a process will be discussed in the next section (it is interesting that a similar 
function was already proposed by Ref. [[FTP). 



On the other hand, it is possible to consider the nonforward distribution amplitude 
( |2-10| ) as an interpolating function between the quark distribution function ( |2.9|) (for Pi = 
P 2 ) and the meson wave function ( |2.4| ), i.e. in the limit P 2 —>■ too. Of course, the 
state obtained by P 2 — > is not the vacuum state, but for mathematical considerations 
concerning the connection of evolution equations for forward and nonforward processes this 
is very useful. As result in section V we obtain relations between different evolution kernels 
(especially between the Brodsky-Lepage kernel and the Altarelli-Parisi kernel) which are 
highly nontrivial. 

Therefore, we can perform two limits 

$ a (x = (1 + t)/2, Q 2 ) = lim q a (t, r, Q 2 ) meson wave function (formally!), 

r— > — 1 

q a (z,Q 2 ) = limq a (z,r,Q 2 ) quark distribution function (formally!). 

(2.13) 

In short we add some further properties of the distribution amplitude ( |2.10| ): 
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Because of the conjugation properties of the light-ray operators (O a (K_;n)) + = 
O a (—n^;n) it satisfies 

( q a (t,T,n 2 )Y = q a (t,-T,v 2 ). (2.14) 
This means that this amplitude is real for P\ = P 2 . 

The normalization of this distribution amplitude follows from the definition ( |2.10| ) and 



J dtq a (t,T,fl 2 ) = {nP + Y l fl V < P 2 \J^(0)\Pl > \nP_=rnP + , (2.15) 

where J"(0) = :ip(0)'~y u \ a il ! {0) '■ is a current with the flavour content a. 

At last we give a representation of the generalized distribution amplitude with the help 
of a "spectral" function which is helpful for later calculations. Taking into account the ideas 
of the general Jost-Lehmann representation for T-amplitudes of two local operators (our 
operators Q2.5|) are in fact bilocal if we use the light-cone gauge in QCD) then the spectral 
functions have a finite support with respect to the new distribution variables z + and z_. In 
this way the matrix elements of light-ray operators can be expressed by " spectral" functions 
f a (z + , Z-, jj 2 ) as follows 



< P 2 \O a {K-;n)\P 1 >= I I dz + dz_e- tK - {hP - )z — iK -^ p+)z+ f a (z+,z.^ 2 ). 

Tli _i_ J J 



(2.16) 



If we insert now this equation ( [2.16|) into the definition of the distribution function fl2.10|) , 
we obtain a new representation 

q°(t, r, Q 2 ) = J J dz + dz. J dK -\™ P +\ e ™-(*P + )(t-z + -Tz-) f a( z+ ^ Q 2 ) 

= Jdz„f a (z + = t-Tz„,z_,Q 2 ), (2.17) 

which shows in which way the mathematically independent distribution variables z + and 
turn over to the more physical variables t and r. 



III. STRUCTURE FUNCTIONS FOR TWO-PHOTON PROCESSES 

One of the first crucial test of perturbative QCD was the investigation of the deep inelastic 
electron proton scattering in the Bjorken region. Theoretically it implies the knowledge of 
the imaginary part of the virtual Compton scattering amplitude. The first QCD treatments 
of this process relied on the light-cone operator product expansion and its relation to the 
moments of the structure functions 



Another important class of processes are exclusive scattering processes of hadrons at large 
momentum transfer 0. Perturbative QCD calculations demand the introduction of wave 
functions describing the incoming and outgoing hadrons. The perturbation theory itself is 
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applicable to a hard scattering part which connects these wave functions. The justification 
of perturbation theory is related to the presence of large momenta. Whereas in the case of 
deep inelastic scattering the scattering amplitude is dominated by the contributions from 
the light-cone singularities in the x-space, here it must not be the case. Therefore its 
perturbative treatment in QCD can be quite different and of course the kinematical regions 
of both processes are also different. 

Here we want to investigate the complete Compton scattering amplitude in the Bjorken 
region near forward scattering fl~9|| . Note that usually perturbative considerations of ex- 
clusive Compton scattering amplitudes are performed in the fixed angle region - were ex- 
perimental data are present or expected [[2(J. If we consider the Compton amplitude in a 
generalized Bjorken region then experimental investigations with present possibilities are 
very hard to realize. On the other hand from a theoretical point of view it embeds the 
Compton amplitude for forward scattering and is therefore interesting in itself. Because of 
the smallness of the momentum transfer between both hadrons it is not possible to introduce 
separately nonperturbative amplitudes (or wave functions) for each hadron. Surprisingly the 
treatment of this process leads in a very natural way to the introduction of new distribution 
amplitudes which for the limiting case of forward scattering turn over into the well-known 
parton distribution functions of deep inelastic scattering. Furthermore these new distribu- 
tion amplitudes satisfy new evolution equations which turn over to the Altarelli-Parisi 
equation for the case of forward scattering. 

To show the virtues of this construction we study additionally the two meson production 
by two photons |[21|| . For small momentum transfer between both mesons also here we 



introduce one distribution amplitude describing both mesons. 



A. Exclusive Virtual Compton Scattering 

Here we consider the virtual nonforward Compton scattering amplitude of a spinless 
particle near forward direction but at a large off-shell momentum of the the incoming photon: 

T^(P + ,P_,g) = i J d'xe^ < P 2 \T (j„ (|) J v (^)) \P t >, (3.1) 

where J^x) = (1/2) : ^(x) 7m (A 3 - X s /V3)^j(x) : is the electromagnetic current of the 
hadrons (for flavour 577(3)). 

The considered kinematics defines a generalized Bjorken region where, in close analogy 
to deep inelastic scattering, this process is dominated by contributions from the light-cone. 

7 *(gi) + H(P 1 )= 7 *(g 2 ) + H(P 2 ), (3.2) 
P + = P 1 + P 2 = (2P,0), E = ^/V + p 2 ), 

P_=P 2 -P 1 = q 1 -q 2 = (0, -2p), q = (l/2)( 9l + q 2 ). (3.3) 

The last expressions inside the brackets denote the values of the momenta in the Breit frame, 
respectively. The generalized Bjorken region is given by 

u = P +q = 2Eq -> oo, Q 2 = -q 2 -> oo, (3.4) 



S 



with the scaling variables 

{ =Z?L, „ = ^ = £^1. (3.5) 
P + q P + q 2v 

To understand this process better we introduce the angle between the vectors p and q in 
the Breit frame by cos0 = -^p-. In terms of these variables we get 

^-(5- H °°^)"> 92 2 ~-(e+^cos0ji/, lywMoos^ (3.6) 

In opposite to deep inelastic scattering the variable £ is not restricted to < £ < 1; for 
example q 2 = demands £ ~ —(\p\/E) cos (p. 

In Appendix A we show that in this region the helicity amplitudes T(X', A) = 
^(A^T^e^A) asymptotically are given by 

A) = (1/2K(A0MA)V (3.7) 

for the transverse helicities, and vanish otherwise. Therefore, only the trace of the scattering 
amplitude has to be considered. 

The first treatment of the Compton scattering amplitude with the help of a (local) light- 
cone operator product expansion OJT8j, symbolized by 



J(x)J(0) = £C„(x a )x* ..2^0^(0), 



n=0 



has the serious disadvantage that the amplitude has to be reconstructed from an infinite sum. 
This may be accepted for deep inelastic scattering where the famous connection between 
the moments of the structure functions and the expectation values of the local light-cone 
operators exists, but in general it is unsatisfactory. The application of the nonlocal light- 
cone expansion p||6|, |T2| , ^6 ,p2| is a suitable possibility to overcome this drawback. In our 



special case it reads (in leading order): 

Jil (I) Jil ~ / dK+dK ~ Fa ^ 2 ' K+ ' K -' ^ ° a ( K +> ^(p 2 ) ( 3 - 8 ) 

with the light-ray operators given in Q2.3|) . The light-like vector 



-/ \ X P 
x{x, p) = x + p- 



x 2 p 2 



p 2 \\ (xp) 2 

is determined by x and parameterized by a fixed constant vector p. The singular coefficient 
functions F a are determined perturbatively; in the Born approximation they are given by 

F a (x 2 , k+, «_) « ie a (2tt 2 (x 2 - ie) 2 ) ^ 5{k+) (<5(k_ - 1/2) - <J(k_ + 1/2)) , (3.9) 

with e a = (2/9)5 O 0+(l/6)5 O 3+(l/6\/3)5 O 8. Possible contact terms arising at the point x = 
lead to trivial contributions to the real part of the scattering amplitude and are suppressed. 
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Inserting the light-cone expansion ( |3.8|) into the trace of the scattering amplitude ( [3.1|) 
and using the representation ( |2.16| ) for the matrix elements < P2\O a (n^\ x)\P\ >, 

< P 2 \O a {K_-i)\P l >= xP+J J dz + dz_e- iK -^ p - )z -- iK - {iP+)z+ f a {z + ,z^P i P ] ^ 2 ) 1 

we get 

T/(P + , P_, q) « 2 j J dz + dz_F a (£, m z+, z.)f a (z + , z_- P t P v fi 2 = Q 2 ), (3.10) 

with 

F a = -e a ftl [ dK _ e i^-(P +Z++ P-^-h-] P + x {6{k _ _ 1/2 ) _ 6 ( K _ + 1/2)) . 
J 2tt z J [x z — ie) z 

(Here the approximation x ~ x has been used.) 

To obtain a useful expression in momentum space the Fourier transform has to be carried 
out. Using 

c ix\q-k\ P + X = ~ k ) _ P +(l k = (p +z+ _|_ P_ Z _) K _ 

2tt 2 (x 2 — ie) 2 (q — k) 2 + ie q 2 — 2qk + ie ' 

(3.11) 

we find 



(3.12) 



In this form the result is unsatisfactory. It depends explicitely on the two distribution 
variables z + , z_ of the distribution amplitude f a (z + , zJ) . However, it is possible to simplify 
this expression further, and to derive an expression for the T-amplitude containing only one 
integral over the standard distribution amplitude. For this purpose we substitute z + = 
t — 7]z_ so that 

T/ « 2 J dt (j^- t - -^—) J dz„e a f a (z + =t- V z-,z-; P % P v f = Q 2 ) 

2 / dt (tt^ - -T-l) e< W, m PiPj\ f = Q 2 ), (3-13) 



where we have used the definition ( |2.17| ) with r replaced by rj, 

q a (t, 77; P t Pj, Q 2 ) = J dz_f a (z + = t - V z„, z_; P^, Q 
Taking into account Eq. ( |3.7| ), we obtain the final result 



T(A', A) « e»(\')e lfl (\) J dt (-L - -L) e a q a (t, m P^-f = Q 2 ). (3.14) 



10 



Note that in order to obtain Eq. (|3.14 ) it was necessary to turn over from the non- 



perturbative function f a (z +} z_) back to the distribution amplitude q a (t } 7]), however, with 
changed variables. The original definition ( 2.10|) of q a (t,r) contains the light-like vector n 



which would be given here by x. This would lead to another definition of the variable r 
by t = xP-/xP + which does not make sense. But instead of this unphysical variable here 
the physical variable rj = qP-/qP + appears. This amplitude generalizes the parton distri- 
bution function of deep inelastic scattering with the above mentioned properties. Note that 
q a (t, Q 2 ) = q a (t, 7] = 0, Q 2 ) is essentially the parton distribution function. 

According to the optical theorem the absorptive part of the scattering amplitude = 
ImT^ I is related to the total cross section in deep inelastic scattering. Since q a (t,r] = 
0, Q 2 ) is real from equation (|3.13|) we get directly 

« - £ dt{6{£ -t)- 6(Z + t))e a q a (t, rj = 0; /i 2 = Q 2 ) 

« " (e a q a (Z, V = 0; /i 2 = Q 2 ) ~ e a q a {-^ rj = 0; /i 2 = Q 2 )) . (3.15) 

Here 

- q a (~^V = 0;/i 2 = Q 2 ) = ?{t,ri = 0;/i 2 = Q 2 ) (3.16) 
represents the antiquark distribution amplitude. 

B. Two Meson Production 

As a further instructive example we discuss the two meson production by two photons 

7*(?i) + 7*fe) = M(P0 + M(P 2 ). (3.17) 

This is a crossed process with respect to the virtual Compton scattering. As variables we 
choose 

q = (l/2)( qi -q 2 ), P + = P 1 + P 2 , P_ = P 2 -P 1 . (3.18) 
Again the generalized Bjorken region is defined by 

v = P + q = 2Eq -> oo, Q 2 = -q 2 -> oo, (3.19) 



with the scaling variables 



{ = Zf, , = £». (3.20) 

P+q P+q 



The virtual photons are taken at large off shell momenta q\ — > — oo, g| — > — oo. Then 
the range of the scaling variable is |£| > 1. In the center of mass system with fixed energy 
of the two mesons these scaling variables can be expressed as 



C-~ q2 - q * + q i ^-Ezl^K^cosS (3 21) 
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In the leading approximation the scattering amplitude can be written as [see ( |3.7| )1 
T(P+, q; A', A) = l -e v 2 {\')e lv {\)i J d'xe^ < P U P 2 \T (j„ (|) r (^)) |0 > . 

(3.22) 

After the application of the operator product expansion, and the calculation of the coefficient 
function we arrive at the same type of expressions as Q3.13j ) 

T/(P + ,P_,g) « 2jdt L±- t - e a q a m (t, m/ j? = Q 2 ), (3.23) 

where the distribution amplitudes reads now 

= / ^Fi^'^ < P^m^xUQ >| ip _^ p+ . (3.24) 



It is interesting that a similar distribution function was already introduced by Chase [17 



in 



the two photon production of two jets. In the limit r\ — > 1 we obtain the known result for 
the production of a scalar meson by two photons for large Q 2 P,|23[|. 



Note that also here we generalize the meson wave function to a two-meson wave function 
in a straightforward way. Of course we cannot predict the dependence of this wave function 
on its internal parameters (this includes also its P1P2 dependence). This is in contrast to 
the large momentum transfer calculation where P1P2 — > 00. Here one introduces for each 
meson its own wave function and then it is clear that the P1P2 dependence of the process is 
included into the hard scattering part. 



C. Parton interpretation 

Here we try to interpret the foregoing results. 

There are two remarkable differences of the foregoing processes in comparison with the 
usually discussed forward or fixed angle scattering processes: 

1. There appear generalized distribution amplitudes describing both hadrons simultane- 
ously. These distribution functions depend, besides other variables, also on the scaling 
variable r\ defined above. 

2. The leading contribution of the hard scattering part consists of an interaction of both 
photons with one and the same internal quark line. The momenta of the other quarks 
are not forced to be large, so that they are confined and can be described implicitely 
by the generalized distribution amplitude. Fig. [I] illustrates this remark for the case of 
the virtual Compton scattering. A standard treatment of exclusive pion (or nucleon) 
Compton scattering at fixed angles would include much more diagrams. 
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FIG. 1. Leading contribution for the virtual Compton scattering in the Bjorken region. 

It remains an interpretation of t in the parton language as momentum fractions of the 
external hadrons. In the case of the exclusive virtual Compton scattering we obtain 

XiP\ + X0P2 

t = (3.25) 

in an infinite momentum frame with respect to P\ + P2- This variable is a natural gener- 
alization of the earlier introduced distribution parameters z or x of the parton distribution 
function and the meson wave function, respectively. We underline that the standard calcu- 
lation of the pion Compton effect at fixed angles for example uses diagrams of Fig. [| 




FIG. 2. Leading contribution for the virtual Compton amplitude at fixed angles. 

Here in opposite to the foregoing case the photons act on both internal lines. Because of the 
large momentum transfer between both hadrons it is possible to introduce separate wave 
functions for each hadron. So the P1P2 dependence is a natural part of the hard scattering 
subprocess. If in addition we can go over to the Bjorken region, then both formalisms should 
be applicable and should lead to the same result. 



IV. DERIVATION OF THE GENERAL EVOLUTION EQUATION 

In this section, we derive an evolution equation for the distribution amplitude ( p. 10 ) 
which contains as limiting cases the Brodsky-Lepage and the Altarelli-Parisi equations. This 
evolution equation exploits the renormalization group equation of the light-ray operators 
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Q2.5|) . For this reason we have to know the renormalization properties of these operators and 
especially their anomalous dimensions. Afterwards we can apply these results to expectation 
values of these operators which are contained in the distribution amplitude. 

It is now well-known that bilocal operators where the distance between both positions 
of the involved field operators is light-like possess additional divergences. Therefore they 
need a special renormalization procedure. In the scalar case the renormalization of light-ray 
operators has been studied from a general point of view by S. A. Anikin and O. I. Zavialov 



12 1 . These operators need special Z-factors depending on the parameters of these operators 
and additional anomalous dimensions [|24[| . A general thorough proof for the renormalization 
of these operators in QCD with covariant gauge fixing we do not know. However practical 
calculations of the Z-factors respectively the corresponding anomalous dimensions exist at 
one loop [|J and two loops for restricted variable ranges [§,[nj for covariant gauge fixing as 



well as in the light-cone gauge. 

It turns out that for a straightforward mathematical treatment the set of operators 
is not well suited. We use, therefore, the more general operators (|2.3|) 



O a (n + , K-]n) = :^(Kin)(n'j)X a U(Kin, K2 , n)'ijj(K2n) :, k± — (/«2 ± Ki)/2. (4-1) 

These operators differ from the original ones by a translation. Their renormalization group 
equation reads formally 



d r 
fi—O a (n + , n) (M 2) = / dK , + dK , _['j( y K + , k' + , k'_; a s (/i 2 )) - 2^(a s (fj 2 )) 



d 

U u {k + , K_;n) (M 2) = / dK' + i 

x5(k + - k' + )5{k + - <)]O a «, n)^. (4.2) 



d\x 

This integral equation contains two integration variables. Using symmetry properties a more 



simple form of this equation will be derived in Appendix B JL6], namely 

, d 
dfi 



d [ 
fi—O a (K + ,K^;n) {ll 2)= / dw + dw-[y(w + ,W-;a s (iJ 2 )) - 27 V ,(a s (/z 2 )) 



x5(w + )5(l — w_)]O a (w + K_ + K + , W^K^.] n)^2\. 

(4.3) 

In the Appendix B we also investigate the anomalous dimensions 7(10+, W-) of the above 
defined operators. As result, we obtain the following support restriction ( |B25| ) 

7(10+, wJ) 7^ for \w±\ < 1, \w + ± W-\ < 1. (4.4) 

Additionally, from the transformation properties of the operators under charge conjugation 
it follows 7(u> + ,u'_) = 7(— w + ,w_) [Eq. (|B32[) 1- 



Technically we investigate in Appendix B the divergent part of all one-particle- irreducible 
(1PI) diagrams with the insertion of one light-ray operator for scalar field theory and QCD 
in light-cone gauge. As essential tool we use the a-representation. However we do not prove 
or carry out the renormalization procedure of the considered operator. For the investigation 
of the support properties of the Z-factors and the anomalous dimensions of these operators 
this is not necessary. So we study the support properties of the divergent term with respect 
to the parameters k + and k_ of the inserted operator and the parameters k' + and k'_ of 
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the appearing counter terms. This allows us to determine the support properties of the 
anomalous dimension of the light-ray operators. 

After these preliminaries we turn to the derivation of the evolution equation. For this 
purpose we differentiate the general distribution amplitude q a (t,r,fi 2 ) with respect to the 
renormalization parameter fi. Thereby we take into account its representation fl2.10| ) in 



terms of matrix elements of the light-ray operator. The differentiation of this operator can 
be performed with the help of its renormalization group equation ( }4.3| ). A straightforward 
calculation runs as follows: 



XS(W+)5(1 — W-)] < P 2 \O a (K-W + , n)\Pi > \hP-=r(hP+). 



(4.5) 



As next step we exploit the translation invariance of the matrix element 

< P 2 \O a {^w+,^w.;n)\P 1 >=< P 2 \O a (K + = 0,K_w_;n)|Pi > e ™-™+(™ p -) (4.6) 
and introduce the new variables k'_ = and t' = (t + tw+)/w- into Eq. (4.5): 



d\x J \wJ\ J 27r(raP 4 

x6{w+)8(l - w-)] < P 2 |O a « = 0,K'_;n)\Pi > \nP_=r{nP + ) 

dt' f (V - t±I?!!±\ [ 7 ( w+) w _) - 2^5(w + )5(l - 

J \W-\ \ W- J 

xq a (t',r,f, 2 ). (4.7) 

As result, we obtain the evolution equation 



(4.8) 



with the evolution kernel 

7(i,t') = J dw-j(w+ =t'w- -t,w-). (4.9) 
From the renormalization group invariance of / dtq a (t,T, fi 2 ) (see Eq. ( |2.15| )) it follows 

dtj(t,t') = 2^, (4.10) 
so that the standard "+" -definition for the generalized function 

[ 7 (t, 01+ = 7(«, - 8{t - / dt"i(t", = 7& - 27**(* " (4-11) 
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arises in a very natural way. We see that instead of the anomalous dimension there appears 
an evolution kernel which contains the anomalous dimension as an essential input. As it 
becomes clear later on it seems to be natural to speak of our general evolution kernel as 
of an extended BL-kernel. The reason for this is the following: Restricting j(t,1f) to the 
parameter region \t\, \t'\ < 1, it coincides with the evolution kernel for the hadron wave 
function, but Eq. (|4.8| ) provides us with a meaningfull definition also outside of this region. 

In this way, we obtain evolution equations for all forward and nonforward matrix ele- 
ments. To our knowledge, in the literature up to now only two examples have been investi- 
gated: the case of forward scattering and the case of the meson wave function. Here, both 
cases are contained as limits P 2 — > P\ or P 2 (see section II). We write these limits down 
in a condensed form: 

• Evolution equation for the quark distribution function (forward scattering): 

q a (z,Q 2 ) = \imq a (z, r,Q 2 ), 

T—>0 

Q 2 ^- 2 q a (z,Q 2 ) = £0jp(^a s (Q 2 ))^(^,g 2 ), (4.12) 



P - 



lim 



-0 2r 



I z z 
7 I r' r . 



Altarelli-Parisi- (Lipatov) kernel. 



(4.13) 



Evolution equation for hadron wave functions: 

This case can be obtained as a formal limit too [compare the remarks after Eq. ( |2.12| )1. 



lim <f(t,r,Q 2 



$ a (x = (l + t)/2,Q 2 ) 

Q 2 ^ a (x,Q 2 ) = f Q dyV BL {x iy] a s {Q 2 ))^{y,Q' 



(4.14) 



Vbl(x,v) = [7(22; — 1,2?/ — l)] + \o<x, y <i, Brodsky-Lepage kernel. 



(4.15) 



Here, the variables x = (1 + £)/2, y = (1 + f)/2 are restricted to < x, y < 1. 



Note that both kernels ( |4. and ( [4.1 5| ) are calculated separately in QCD. From the above 
considerations it is obvious, however, that they have a common origin, namely the anoma- 
lous dimension r y(w + ,W-) of the light-ray operators being hidden within the corresponding 
amplitudes. We return to this point in section V. 

Let us add a remark concerning the distribution amplitude introduced for the nonfor- 



ward virtual Compton scattering. This distribution amplitude appearing in ( 3.14 ) or ( 3.23 ) 
satisfies in the flavour nonsinglet case the evolution equation 



1 dt' 

-iSI 




(4.16) 
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Here, the distribution function depends on variables with a clear physical interpretation. It 
coincides formally with the definition (4.S) so that no further proof is needed. Nevertheless, 
this equation can be proved directly: It is possible to derive an evolution equation for 
f a (z + , z_,Q 2 ) starting from the original renormalization group equation, and to use Eq. 
( p. 17] ) for the amplitude q a (t', rj, Q 2 ). This chain of relations is independent of the concrete 
physical meaning of the variables, so that the evolution equation (|4.8|) is valid for the changed 
variables too. In the singlet case we have a mixing problem with a gluonic contribution. 
This problem could be solved in a straightforward manner, too. 



V. THE EXTENDED BL-KERNEL, RELATIONS 
BETWEEN THE BL- AND AP-KERNELS 

In section IV we derived the general evolution kernel (an extended BL-kernel) for distri- 
bution amplitudes containing the (restricted) BL-kernel and the AP-kernel as special cases. 
In QCD the BL-kernel is calculated as an evolution kernel for a restricted range of the vari- 
ables. The interesting questions are: Exists there a possibility to determine the extended 
BL-kernel knowing this restricted kernel only? Furthermore, if such a procedure exists, is 
it then possible to apply it to the existing two-loop QCD calculations? If yes, then one 
should be able to perform the limit to the Altarelli-Parisi kernel which should then coincide 
with the already existing results obtained by straightforward computation. We solve these 
problems in the following three subsections. 



A. General Aspects 

Here, we discuss the general solution of the problems stated above. Starting from the 
support properties of the anomalous dimension 7(u>+, and from the definition of the 



evolution kernel ( O ), we study its domain of definition. Afterwards, we study the exten- 
sion procedure. It turns out that the restricted BL-kernel contains already the essential 
information and that an explicit continuation procedure can be prescribed. 

The first question is: What is the correct region in the (t, t')-plane where the kernel 
7(t, t') is defined in fact. To answer it, we start from the general representation ( |4.9| ) of this 
kernel 



j(t,t') — J J dw-dw + 5{w + — t + t'w-)j(w + , W-) 

and use all known symmetry properties and support restrictions: 

j(w + , W-) = 7(— w + , W-), \w±\<l, \w + ±W-\<l. 

After some algebra we obtain the following representation in the (t, i')-plane 0] 

7 (t, = [6(t - 00(1 -t)~ 0(f - t)9(t - l)]f(t, 

+[9{t' - t)6(l +t)- 6(t - t')9{-t - 1)]/H, -0 
+[0(-t - 00(1 +t)~ 0(t' + t)0(-t - l)]g(-t, 
+[9(t + t')9(l -t)- 9{-if - t)9(t - l)]g(t, -0, (5.1) 
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where the functions f(t,t') and g(t,t') are given by 

f(t,t')= 1 t dw^j(w + — t — W-t' , 



g(-t,t') 



i+t 

l-i' 



<iw_7(w_| 



— t — W-t , — iu_ 



(5.2) 




FIG. 3. Support of j(t,t'), where f±± = f(±t,±t'), and g±^ = g(±i, =Fi') are defined by Eqs. 
P). 



Note, that for \t\, \ t'\ > 1 because of support restrictions the lower boundary in the integral 



representations ( |5.2| ) is not attained. However, in these regions only the differences 

f(t,t')-f(-t,-t') and g (-t,t')-g(t,-t') 



appear in Eq. (|5.1|) , so that such undetermined contributions eliminate each other. It is 
clear, therefore, that the evolution kernel is defined in the complete (t, £')-plane as shown in 
Fig. |. 

The second question is now: Assuming we are able to calculate the anomalous dimension 
in the region —1 < t,t' < 1, is it possible to determine it in the whole (t, t')-plane? This 
problem contains the principal question whether the restricted BL-kernel contains already 
the full information on the general anomalous dimension. For solving this problem we 
perform a Fourier transform of the extended BL-kernel and later on of the restricted BL- 
kernel. It will be shown that because of the holomorphic properties of these functions it is 
possible to identify them. 
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So we perform a Fourier transform of the extended BL-kernel. Thereby, we take into 
account the known support restrictions explicitly 



7(A,f)= /<fce" t [ 7 (t > 0]+ 

r dw + dw^{w + , w-) (e lX{w+ +t ' w ' ] - e iXt ') . (5.3) 



I W\ | , | w_ | < 1 

I W<_ ±KJ_ I <1 



Because of the finite integration region in the Fourier integral with respect to the variables 
w + , w_ the resulting Fourier transform of the anomalous dimension is an entire function of 



the variables A and t' 25 



Now, we define a second entire function of the same variables by the Fourier transform 
of the restricted BL-kernel 

V(\,t') = J dte iXt [V{t,t'))+ 

dw + dw^9(l - u>_t' - w + )0(w-t? + w + + 1) 

x 7 (w+, wJ) [e l ^ w + +t ' w -^ - e iXt> ) . (5.4) 



1 | , \ w_ | < 1 

\w_l_ I < 1 



Both functions have to be compared. The essential point is that, due to the known support 
properties, for \t'\ < 1 it follows automatically also \t\ < 1 so that for \t'\ < 1 there exists 
a region where both expression j(\,t') and V(X,t r ) are identical. Therefore, knowing the 
Fourier transform of the BL-kernel we can obtain the Fourier transform of the extended 
BL-kernel by analytic continuation. After performing the reversed transformation 

r d\ 

[7(M')]+= / ir-e- iXt {V(X,t')} AC (5.5) 

we obtain j(t,t') completely. {V(X,t')}Ac denotes formally the analytic continuation pro- 
cedure to be carried out. 

Consequently, it is sufficient to perform all diagrammatic calculations in the restricted re- 
gion \t\, \t'\ < 1. The essential information on the anomalous dimension is already contained 
inside of this region. 



B. The Extended BL-Kernel of QCD in Two-Loop Approximation 

As an exercise to the foregoing investigations we will extend the Brodsky-Lepage kernel of 
QCD into the complete (t, t')-plane starting from the explicit one- and two-loop calculations 
of the restricted BL-kernel PJlQ|. 

The known calculations of the evolution kernel P,|T0f are performed in the range — 1 < 
t,t' < 1, which was sufficient for studying the evolution of the hadron wave function. The 
difficult two-loop calculation have been done in the variables x — (1 + 1)/2, y — (1 + 1')/2. 
So for the rest of the paper we shall prefer these variables too. The result of the above 
mentioned calculations is 
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V(x,y) = ^c F [V (x,y)} + + (^J c F [c F V F (x, y) + (c G /2)V G (x,y) 

+ (N F /2)V N (x,y)] + + 0(c 



with 



V (x, y) = 6(y - x)F(x, y) + 



V N (x, y) = 0{y - x)v N (x, y) + 



x — > x y 

y^yf 

j x — > X 



Vg(x, y) = 0(y - x)u G (a:, y) + y) + 



y^y . 



Vf(x, y) = 6{y - x)v F (x, y) - —k(x) - G(x, y) + < 

^y [y ~~ ^ y , 

Here, the following auxiliary functions are used: 

F(x,y) = -[l 



y V x-y 



10 2x 2 fx 
v N (x,y) = --F-----Fln f- 



, x 67^ 17 x 11 , fx\ 
v G (x, y ) = — F + — - + — Fin - , 
9 3 y 3 V J// 



3 1/2 \yj \ x-yj \y J 



k(x) = ln(x)(l + ln(x) -21n(x)), 
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G(x, y) = 2d(y - x) [FLi 2 (x) + F ln(y) ki(x) - FU 2 (y)] +6(x + y-l) 



2(F - F) 



xLio [ 1 - - 



(F - F) \n 2 (y) - 2F\n(y) \n(x) + 2FU 2 (y) - 2FU 2 (x) 



where the Spence function Li 2 is defined by 

Li2(x) = _r dz Mi^£). 

Jo z 
As shorthand notation we have introduced 

x = l-x, y = l-y, F = F(x,y) = F(l - x, 1 - y). 



(5.8f) 



(5.9) 



(5.10) 



The group theoretical constants are cp = 4/3, cq = 3. N F denotes the number of quark 
flavors. Contrary to the usual notation of the kernel in Ref. [pl,[l0|| we have changed the 
notation by 

G(x, y) + 2% - x)F\n(y) ln(x) -> G(x, y). 

Furthermore, for convenience we have introduced the function k(x), and in the last row of 
( |5.8d| ) we have used the decomposition F(x,y) = — | + 

The obtained kernel V(x, y) is not defined for x = y, however the "+" -definition provides 
a regularization prescription. 

The kernel ( ^.6|) has to be extended using the foregoing results. Let us begin with the 
results of the one- loop calculations ( p.7a| ) and ( |5.8a| ). According to our procedure, we have 
to perform an analytical continuation of the Fourier transform 



V {\y) = 2 
which may be written as 



dx 



6(y - x)F(x,y) + 



\ x - 


-> X 




\y- 







(5.11) 



V (X,y) = 2 f y dxF{x,y) (e^ 2 *" 1 ) - e'^ 2 ^ 1 ') 



+ 2 dxF(x,y) (e 



3 iA(l-2x) _ JX(1- 



(5.12) 



Both integrals are well defined and represent an analytical function with respect to the 
variables A and y. In the first integral we introduce the variable z = x/y and in the second 
the variable z = x/y so that 



V (X,y) = 2 f dzyF(zy,y) (e^ 2 ^ - e iX ^) + 



(5.13) 



Obviously, yF(yz,y) = y — 1/(1 — z) is an entire function with respect to y, and the 
integrand yF(zy,y) (e^ 22 ^ -1 ** — e l ^ 2y ~ 1 '\ has no singularities for < z < 1. Then the 
analytical continuation of Vq(A, y) is given by 
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{%{\v)}ac = 2 jf 1 dz{yF(zy,y)} AC (e^ 2 ^ - e iX ^) + {*^\} , (5-14) 

where {yF(zy, y)}AC — V ~ 1/(1 — -2) denotes the analytical continuation with respect to 
the variable y. Obviously the inverse transformation provides 



^(x, y) = /' ^{yF(zy, y))}^^ - 2zy) - - y)) + 
Jo 

= e(l-z)e(z)sign(y){F(zy,y)} AC \ z =z 

-5(x-y) dz{yF(zy,y)} AC 
Jo 



x — > X 



rr — > rr 
. 3/ — »- 27 . 



6(1 - z)0(z)siga(y){F(zy,y)} A c\z=z] + + |* ^ * J . (5.15) 

So we obtain as extended BL-kernel in one-loop approximation 

V„"'(i,!/) = 7o(« = 2i - 1,C = 2j/ - 1) 

= S (l-£) () (£) s i g n( !/ )£(l-^ N ) + {^!}. (5.16) 

V yj \yj y\ x-yj yy^y) 

Note that this result can be easily checked by a direct calculation of 70 (i, t'). 

The same procedure can be performed for the two-loop contributions to the BL-kernel. 
The calculations are much more complicated. For convenience, we consider each term sep- 
arately. 

We start with the contributions coming from the functions , v G and v F . They 
have - as in the foregoing case - the property that the functions yvN(zy,y), yvc(zy,y) and 
yvp(zy,y) are entire analytical functions of the variable y. Therefore, the extension of the 
corresponding part of the BL-kernel can be performed in the same manner: The extension 
of these contributions is obtained by substituting 

% - x) - 0(l - Q sign(y) (5.17) 

and performing the analytical continuation of the functions vjy, vq and vp, e.g. 

v F (x,y) = {v F (zy,y)} A c\*=*- (5-18) 

More problems occur for the contributions J=M§) an d ^jk(x). These functions have a 
pole for y — 1. However, the Fourier transformation of 

«■(*, ») = •&-- (I) - + (5.19) 

can be written in the following manner 

K(\, y) = f G dzUy 2 k(z) (e^- 1 ) - e^" 1 )) - {y = 1}] + { } . (5.20) 
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Obviously, the integrand in Eq. ( |5.20| ) is an entire function, and we can use the procedure 
described above. 

The extension of G(x, y) is more complicated. If we write down its Fourier transform 
in the same manner as in Eq. (|5.13| ) then the integrand seems to contain poles and cuts. 



These are caused by the F-function and by the In- and Li2-functions, respectively. However, 
with the help of identities for the Li2-functions it is possible to find such a representation 
in the region < y < 1 that the integrand is an analytic function (without poles and cuts) 
for < Re y. Similarly, it is possible to find a representation where the integrand is an 
analytical function for Re y < 1. With other words, the Fourier transform of G(x,y) is 
in fact an entire function, but, opposite to the cases treated above, there does not exist a 
representation for the integrand in terms of Li2-functions which is well suited for all values of 
y. Therefore, formally we can obtain the extension by performing the substitution ( p,17| ) (of 
course, 6(y — x) — > 9(1 — |)#(|)sign(?/)) and the analytical continuation of the corresponding 
functions, 

G ext (x,y) = 2e(l-?)e Q sign(y){FLi 2 (x) + F\n(y)\n(x) - FLi 2 (y)} AC 
+6 U-fjO (~) sign(y) |2(F - F)Li 2 (l - ~ 





-2Fln(y) ln(x) + 2FLi 2 ( 

As result, we find for y < and y > 1 terms which contain an imaginary part. However, 
these cancel each other exactly so that in the physically interesting regions the extension 
provides in fact a real function. 

As final result the extended BL-kernel reads: 

[ 1 (2x-l,2y-l)] + = V ext (x,y) 



^-c F [V ext (x,y)] + + c F [c F V^{x,y) + (c G /2)V^\x,y) 

+ (N F /2)Vff :t (x,y)} + + 0(a 3 s ), (5.22) 

where the extension of Vq, V f Vq and Vn (|5.7a| - |5.8^) follows from the discussed procedure 



which is given by ( |5.17| ), ( |5.18| ) and ( p.21| ). So, up to the changed ^-structures the analytic 



expressions are formally identical to the originally calculated ones. 



C. Consistency between the QCD Calculation of the BL- and AP-Kernels 

An interesting application of our formula ( 4.13|) would be a check of the mutual consis- 
tency of the very complicated two-loop calculations of the BL-kernel and the AP-kernel 
||. As it has been shown the extended BL-kernel should contain the AP-kernel as limiting 
case. Therefore, independent calculations of both kernels are not necessary. Otherwise, if 
they exist, then they must be consistent with each other. 

The existence of a connection between the AP- and the BL-kernels follows already from 
a very simple argument. Indeed, looking at the local anomalous dimensions it seems to be 
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obvious that the general anomalous dimension matrix V nm = 7„ m — 2^5 nm corresponding 
to the BL-kernel \fM] 



f dxx n V BL (x, y)=jz V nm y m (5.23) 



„1 n 

I 

m=0 

contains the diagonal matrix elements corresponding to the AP-kernel, 

dzz n P(z) = V nn , P{z) = — / dnz- n - l V nn . (5.24) 

If the local anomalous dimensions would determine the corresponding kernels then it should 
be possible to determine the AP-kernel from the given BL-kernel. Starting from this obser- 
vation there were several incomplete trials to solve this problem [[J. 

Here, we shall determine the AP-kernel from the known extended BL-kernel. In the 
language of the local anomalous dimension we can say that in this way an independent 
check for the correctness of the diagonal anomalous dimension will be performed. A check 
for the nondiagonal elements of the anomalous dimension matrix has been carried out by 
exploiting consequently conformal symmetry breaking [[11]]. Therefore, the considerations 
given in the following enable us to conclude that the two-loop approximation of the BL- 
kernel has been correctly computed. 

To obtain the AP-kernel practically, we have to apply Eq. ( |4.13| ) which we reproduce 
here for convenience 

P{z) = lim M- 1 ^ - 1,23/- l)]+\x=**=$ (5.25) 

r — >U 

Let us remark that according to this equation it is just the added new region (|t|, \t'\ > 1) 
which is essential for the determination of the AP-kernel. 

For technical reasons, we first discuss the ^-structure and then we consider the limiting 
procedure for the coefficient functions contained in Eqs. ( |5.8aj - ^8^) . The limits r — > with 
x = z/t, y = 1/t for the ^-functions are 



0(^)0(1-^ 

0(^- y )0(l-i- y 



(z)9(l-z), (5.26) 
{-z)6(l + z). (5.27) 



It is important for further calculations that "related ^-structures", i.e. structures which turn 
into each other under x <-> x = 1 — x and y y = 1 — y - have the same limit. We note, 
that all expressions have the following typical structures 



S 1 (x,y) = B\^Je(l- BigaWTfay) + 9 (J J 9 1 1 - |J sign^x, y) 



(5.28) 



or 
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S 2 (x,y) = $( ^je (l - sign(y)T 2 (s,y) + f^j f 1 - \ sign(y)T 2 (x,y). 



(5.29) 



So we have to determine the following limits: 

limlrl- 1 ^ (^i)=^(l-,)limr- 1 {T 1 



z 1 
r r 



z 1 
r r 



(5.30) 



and 



lim Irl 1 So f -, - 



r t 



.,)«(! + ,) tar- {T,(£,i)-T,(l -£,!-!)} 



As short notation we introduce the symbol LIM, 

LIMT(x,,)^hmr-{T(^I)-T(l-^l-I)}. 



(5.31) 



(5.32) 



The determination of these limits will be performed in Appendix E. 

Putting together all results obtained there we get an expression for the AP-kernel fol- 
lowing from the extended BL-kernel ( |5.6| - |5.81{ ). In analogy with Eq. ( |5.6| ) we write 



P(z) = ^- CF [P (z)] + + c F [c F P F (z) + (c G /2)P G (z) + (N F /2)P N (z)}. 

+ OK 3 ). 

Let us compute each term separately 

P Q (z) = ]m k \T\- 1 V (-,- 



(5.33) 



[z)d{l-z) 



1+z 2 
1-z 



(5.34) 



according to ( Elb|) . 



P N (z) = \im\r\- 1 V N (-,- 



(z)0(l-z): 



l±f!(.„ W + !) + 2(l-,) 



(5.35) 



because of (|Ela| - [Elc 



P G (z) = ^m\T\- 1 V G (-,- 



(z)6(l - z) 



1 + z 2 / 2 . . 11 . . 67 tt 2N 
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+ 2(l + z)ln(»+-(l-z) 



+ 0(-z)O{l + z) 



1 + \z[ 
' 1 + \z\ 



Li 2 





z 




1 + 


z\ 



Li 5 



1+ z 



+ \^(\z\) 



ln(|z|) ln(l + \z\) + 2(1 + \z\) In(|z|) + 4(1 - \z\) 



(5.36) 



Here we have taken into account ( [Ela| - |Elc]JElg]jETlT]) . For the last term we obtain the result: 

'z 1 



P F {z) = limlrl^Vp 

r— >0 



T T 



(z)9(l - z) 
4 



1 + f \n 2 {z) + ^ ln(z) + 2 ln(z) ln(l - z) 



l-z 



1 4- 3^ 2 

\n 2 {z) --{?> + 7 z)\n{z)- 5(1 -z) 



2(1 -z) 



1 

2 ( 



-0(-*)0(l+*) 



1 + \z\< 

'TTR 



Lio 





z 




1 + 


z\ 



Li 2 



1+ z 



^ln 2 (kl) 



lu(|*|) ln(l + \z\) + 2(1 + \z\) ln(|z|) + 4(1 - \z\) 



(5.37) 



The obtained result for the AP- kernel ( |5.33| - |5.37| ) is completely equivalent to that of Ref . || . 
It also includes in a natural way the second order contributions stemming from the internal 
ant i- quark lines. 

Let us add a remark concerning the regularization prescription of the considered kernel. 
It is very important, that our regularization prescription for the extended BL-kernel turns 
after performing the limiting process, automatically into the well accepted regularization 
prescription for the AP- kernel. If we start with (|4.11|) and ([5.25 ) 



lim It 



-1 



Z 1 

71-,- 



T T 



Z 1\ f°° , , / , 2 

/ dx'q \2x' - 1, - - 1 

T T J J-oo V T 



(5.38) 



then the limit of the ^-structure restricts the integration range in Eq. ( ET3SD . The substitu- 
tion x' = z'/t leads to 



lim |t| 1 



7 (^) l™J z 'l{^\) =P(z)-6(z-l) Jdz'P(z') 



(5.39) 



This completes our considerations. 
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APPENDIX A: VIRTUAL COMPTON SCATTERING AMPLITUDE 
IN LEADING APPROXIMATION 



The aim of this Appendix is to study the behaviour of the helicity amplitudes of the 
Compton scattering amplitude in the generalized Bjorken region. The helicity amplitudes 
are defined by 

r(A',A)=^(A')T^(A), (Al) 

where £,(A) denotes the polarization vector of the virtual photon (z = 1,2): 

e?(A) = (0,ei(A)), r<(A)g< = 0, e i (A)e i (A / ) = 8x\>, for A, A' = 1,2, 

e?(A = 3) = |g?|- 1 / 2 (|5| ) g?g/|gS|) > 

= \li\~ 1/2 (Qi - QiCi/ciqi), cf = (1, 
ef(A = 0) = l^- 1 /^. (A2) 

It will be shown that in leading order they can be expressed as 

i4(A')^(A)T/ for A', A = 1,2, 



v ' ; \0 otherwise, v ; 

so that in section III we need to consider the more simple expression T v v only. 

To get these results at first we remark that because of the current conservation c^T^ = 
and T^q" = the helicity amplitudes T(A', A) vanish for A' = or A = 0. For the same 
reason a kinematical suppression of the longitudinal degrees of freedom occurs, for example 

\4\ 



4(A = 3)7^ = \q 2 2 \- 1/2 q%T, v ± ^^c^ u = ±^1^ for q 2 < q (A4) 
is suppressed because the factor 



to±p-m_.«F (_l_\ i = 1 , 2 . (A5) 



Ciqi Ci(q±P-/2) dq \VW , 

vanishes asymptotically. 

For further simplification we need explicite representations and approximations of the 
scattering amplitude. Because of our kinematical assumptions, in x-space the scattering 
amplitude ( |3.1| ) is dominated by contributions arising from the neighbourhood of the light- 
cone. So it is possible to apply the light-cone expansion of the current product. Instead of the 
standard expansion in terms of local operators we use the more effective expansion based on 



light-ray operators |12| , |22| . This expansion relies on a perturbative expansion of the product 
of the electromagnetic currents, too. The crucial difference between both expansions is the 
following: Whereas usually the light-cone expansion follows from a suitable Taylor expansion 
leading to the "local" light-cone operators, here a corresponding Fourier transform generates 
the nonlocal ("light-ray") operators. In leading order we obtain |4j] 
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/ '-/,, ( % ] J v ( — ^ I PS y dK, + d,K-F a (x 2 , K+, /J 2 )r^X a Op(K + X, K-X)t,,- 



2 



(p 2 ) 

+ / dK + d,K^F a (x 2 , k + , k_, fi 2 )e^^x a O^(K + x, K-x)^, 



(A6) 



where T^ 13 = g^ u g af5 — g^gj 3 — gj^gj 1 and e Q/3 denotes the e-tensor and 

xp 



x(x, p) = x + p- 



x 2 p 2 



P 2 \\ ( X P) 2 

is a light-like vector determined by x and parameterized by a fixed constant vector p. 
In the tree approximation we obtain for the coefficient functions F a ( ) and 

F a (x 2 } K + ,K_,p 2 ): 

F a (x 2 , k + , «_) = - , - 6(k + ) (6(k- - 1/2) - + 1/2)) 
Z7r z (x z — ley 

F (x 2 , k + , «_) = 27r2(x2 _ ze)2 ^ + ) (*(*- - 1/2) + + 1/2)) , 

2 1 1 
c a = -S a o + + w^^s for flavour SU(3). (A7) 
9 6 6V3 

In this approximation the gluon operators are automatically suppressed, and only the quark 
operators contribute: 

Oq(k + x, K-x) = : i})(K\x)^p\ a U{K\X, K2x)i/j(k2x) :, (A8) 
Op(K + x, K-x) = : ip(Kix)i , y r y^X a U(KiX, K2x)ip(n 2 x) :, «± = (« 2 ± «i)/2. (A9) 

Here, U(k\X,K2x) is the standard phase factor ( |2.5| ), and A a determines the flavour con- 
tent, additionally we defined A = 1 so that the summation over the index a includes also 
the flavour singlet case. The choice p 2 = Q 2 for the renormalization point justifies the 
perturbative calculation of the coefficient functions. 

The nonperturbative part of the Compton scattering amplitude is contained in the matrix 
elements of the quark operators. Leaving aside the nonessential /t + -dependent exponential 
[in fact it drops out because of S(k + ) in Eqs. (|A7|) 1 

< P 2 \0%{k+x, >=< P 2 \0%(k-x)\Pi > e iK+( * p - } 

we write down the general structure of the matrix elements of these operators between scalar 
or pseudoscalar mesons or spin averaged baryon states 

< P 2 \TO*(k_x)\P 1 > = P + ^(K_x,P,,/i 2 ) + P^t a 2 (K_x,P l ,p 2 ) 

+ K_^M 2 C(/t_x,P i , y u 2 ), (A10) 

< P 2 \t6%(k..x)\P 1 > = e^SeP+P-K-X e i a (K_S;,P h p 2 ). (All) 

The mass factor M 2 = P 2 in the third term of Eq. ( |A1(J| ) has to be introduced for dimensional 
reasons. This term does not contribute in leading order. Inserting Eqs. (|A6|) , ( |A10|) , ( [All 
into Eq. (|3.1|) we are able to write formally 
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T^(P+, P_, g) = (pP-£-T 1 (P + , P_, g) + At 2 ( P+) g ) 



(A12) 



where 



Ti(P + , P_, g) = /" d 4 x J dK+dK^e ixq+iK+isp - ) F a (x 2 , «+, P, Q 2 ). (A13) 



A corresponding formula is valid for T(P + , P_, g). The r.h.s. of Eq. (|A12 ) can be simplified 
further: Since 



dO 1 



qP + <9£ qP + drj 



(A14) 



we can estimate 




dq c 



T i (P + , P_, g) « -2g" [ + -L *L ) T ;(£, ^ Q 2 ) 



<9Q 2 + gp + 

Q 2 ^+^)m, V ,Q 2 ) for Q 2 ,gP + ,gP_^oc. 



2g c 

Q 2 l«» ' 



(A15) 



Using 



e J (A)g = -^(A)(g J + (-l) i P_) = -(-l)^(A)P„ for A = 1,2, 



(A16) 



and the explicit expression for T^ ua p, we obtain after some algebra the leading contributions 
for the transversal helicity amplitudes: 



T(A', A) « e%(\')e lfl (\) (p+^T, + P_^T 2 N 



(A17) 



Here we took into consideration, that the T amplitude gives asymptotically nonleading 
contribution, which follows from similar arguments. Using Eq. (|A12), a direct calculation 



of the trace shows 



(A18) 



Comparing Eqs. ( |A17|) and ( |A18|) , we obtain the desired result ( |A3| ) . 
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APPENDIX B: PROPERTIES OF ANOMALOUS DIMENSIONS 
OF LIGHT-RAY OPERATORS 



This Appendix is devoted to a general discussion of the support properties of the anoma- 



lous dimensions of light-ray operators (|2.3|) 



O a (K + , n) =:ip{Kin){n^)\ a U{Kin^ K 2 n)^(K2n) :, k± — («2 ± «i)/2, n 2 = 0. 

(Bl) 

These operators differ from the original ones ( |2.5| ) by a translation. Their renormalization 



group equation reads |L6| 
d 



^0\k- fi)^) = / dV ( 7 (fi, ^; t^ 2 )) - 2 7 ^(/i 2 ))<5( 2 ) (k - «0) O a (tl; n)^y 



(B2) 



For convenience we divide the anomalous dimension of the operator into the anomalous di- 
mension of the 1PI vertex function (with two external momenta) ^(k, k/) = 7(k+, k'_) 
and a part which is proportional to the anomalous dimension of the quark field 7^ . As a 
shorthand notation we use O a (n;n) = O a (n + , n_;n), d 2 K/_ = dn' + dK'_ and 5^ 2 \k — = 
<5(k+- k' + )8{k'_- 

In the following, we investigate the anomalous dimensions j(k,is[) of the above defined 
operator using the a-representation for 1-particle-irreducible (1PI) diagrams containing this 
operator insertion. As result, we shall obtain the support restriction ( |B25| ), i. e. \w±\ < 
1, |w + ±w_| < 1 with w + = (k' + — k + )/k-, = k'_/k- and the special variable dependence 
— ^Z 2 l{w +) wJ). Additionally, from the transformation properties of the operators 
( P30| ) under charge conjugation it follows r y(w + ,w_) = r y(—w + ,w_). In all what follows we 



give the more technical proof of these results. 

In principle, we have to investigate the renormalization of gauge invariant operators on 
the light-cone. Thereby we should have in mind: 

1. It is essential that these nonlocal operators are defined on a light-ray. This changes 
their renormalization properties and induces a close correspondence to the local oper- 
ators known from the standard light-cone expansion [B7J. 

2. From the renormalization procedure of gauge invariant operators it is well-known that 
during the renormalization process there may appear gauge variant and ghost oper- 
ators. In the complete renormalization matrix they appear in a triangular form, so 
that the anomalous dimensions of the physical operators remain unchanged. Moreover 
these operators have vanishing matrix elements between physical particle states. For 
this reason, we do not need to consider such contributions 



3. Up to these restrictions the operators ((Bl]) form a complete operator basis of minimal 
twist and minimal dimension. 

The anomalous dimension can be determined from the 1PI vertex function of the operator 



to be studied For instance, by using dimensional regularization it can be shown that 



in the minimal subtraction scheme the following simple rule is valid [29 
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d 

T(«> £ 9) = 9q~ Z 9) (B3) 

where can be obtained from the residue of the 1PI vertex function R'G(pi,p 2 ; g\ £\O a (n)) 
with respect to the parameter e = (4 — n)/2. In detail, by R'G(pi,p 2 ; g\ e|O a (/t)) we denote 
the 1PI vertex function (with two external momenta p\ and p 2 ) containing the operator 
insertion where all subdiagrams are renormalized, and only the overall renormalization is 
not carried out. The e-expansion of this function reads |30j 



R'G(p 1 ,p 2 ;e\O a (K;h)) = G^\ Pl) p 2 \O a (^n)) + -G^{p u p 2 \O a {K;,n)) + -•-. (B4) 

In general, the 1/e term is given as an integral with respect to the k- variables of the Z-factor 
with the vertex 

0^{K;n;p) = (f l7 )AV K+( ^+ )+iK - ( " p - ) (B5) 
of the bare operator (|B1| ). We have, therefore, 

G [1] (p l ,p 2 \O a {! 1 ;n)) = J d 2 tlZ^{K,!l)0^{tl;n;p). (B6) 

Using the special form (|B5|) of the operator vertex we can write 



(B7) 



{fi 1 )\ a Z^{K,tsL)= /'^^e-H(^ + )--'-(%-)G'W(p 1 ,p 2 |O a (K;n)). 
J 2,71 J 2n 



Because of Eqs. ( p3|) and (|B4]) the investigation of the support properties of the anoma- 
lous dimension requires the knowledge of the support properties of all possible diagrams 
contributing to R'G(pi,p 2 ; g; e|O a (V)). 

The standard method for the investigation of Green's functions is the application of the 
a-representation. We use the following formula for the propagator 



lim p [ LI. ) i f°° dae ia(k 2 - m ^ P )^k , B8 n 
£-0 \2id£ Jo ' K ' 



w? — k 2 — ip €^0 y2i d£ 

where P(k) denotes a well spezified polynomial in the momentum k. Also a complicated 
matrix structure of the vertices can be included into this polynomial. With these remarks 
the a-representation for a 1PI diagram with L internal lines and V vertices and two external 



legs can be written as [22 



r°° ( 1 8 \ 

Or(pi,p a ) = j Q ^n^^^J e "" i(mi " ip) /r(a;p;i)li=o, (B9) 



where 



const. i{E' J=1 P l ^(°)P J +2^f =1 ^^ 1 € i B ! aa)p I -Efc i =i^ X w(^} 
M«,MJ - D n/2( a f 



(BIO) 
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Thereby, the regularization of the ultra-violet divergencies is ensured by dimensional reg- 
ularization, i.e. by an analytical continuation from n = 4 to n = 4 — 2e, e > 0. For later 
convenience, we have displayed the polynomial structure with respect to variables Pi and £i 
explicitly and introduced coefficients depending on a-parameters only. Thus, our notation 



differs slightly from that of Ref. |22| ; furthermore, in contrast to Ref. |22] we have absorbed 
the overall factor 1/D of the exponent into the definition of Aij, Bu and K^i- Anyway, the 
coefficients A^j, Bu and K^i may be simply deduced from the general rules given there. 

For our purpose we need diagrams containing one insertion of the light-cone operator 
(PT|). Of course, the vertex corresponding to this operator is complicated. However, if 



we consider QCD in axial gauge n^A^ = 0, then the phase factor U(nih } K2n) does not 
contribute, and only the bare vertex (|B5|) remains. Moreover, ghost operators do not exist. 
The drawback of this procedure is the comparatively complicated structure of the gluon 



propagator ]3T 



—i 

ab = °ab-- 



X = 5, 



k 2 + ip 



, „ n^k u + n^k^ 
9^ . 



nk + ir nk — it 



(Bll) 



The fi-dependent terms contained in the second part of this propagator need additional 
considerations. For simplicity we choose for the regularization of the unphysical pole term 
(nk) -1 the principal value prescription. It is well-known that this procedure generates non- 
local counter terms. This problem can be avoided by using the Mandelstam or Leibbrandt 
prescription p2| . However, in our case this is not necessary, because the considered physical 



quantities are independent of the gauge fixing procedure and therefore, do not depend on 
the prescription of the unphysical poles. As usual, the spurios divergences are subtracted 
by hand (see also the two-loop calculation in Ref. |§). Below we shall show that the con- 
tributions of the unphysical poles do not change the support property of the anomalous 
dimension. So we shall ignore them; in fact this means that we study the scalar theory first. 

The operator vertex ( |B5| ) may be included into each diagram as a special first vertex 
connected through the lines 1 and 2 to the remaining part of the diagram. In the a- 
representation, we have 

G r (p|0(«;n)) = O v ^«;n; G r (p,£)||=o, (B12) 

with 

O v (K;n;p) = e ^+(«p+)+«-(™p-). 

The differential operator in the exponential shifts the variables according to £i — > ^i + ^n and 
£2 — ► £2 + ^n. This shift can be taken into account explicitly, so that the a-representation 
with the special operator insertion takes the form 

roc L 



G T (v\0^n)) = J~ ' daft Pi f-^-j e^^/rfeBS 



1=1 

x exp i < KiBijinpj) K i K n(^i) \ \i=°- ( B13 ) 

I i,j=i »=i 1=1 
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Here we have used n 2 = and the symmetry relation K k i(a) = Ki k (a) which is proved in 
Appendix C. After summation over all allowed graphs we obtain the 1PI vertex function 

G(p\0( S n)) = G t (p\0(k; n)). (B14) 
r 

Let us now turn to the support properties of the anomalous dimension with respect to 
the variables k and k/_. We use the fact that the partly unrenormalized 1PI vertex func- 
tion G(p\0(Kr,n)), the renormalized functions RG(p\0(tk, h)) and the anomalous dimension 
1{!5li!5l) have the same support. So we study the support of the vertex function ( B14|) . Its 
support restrictions follow from the properties of the K-dependent part at £1 = ... = £l = 
in the exponential: 

2 

/V+ [(fin + Baa + B 2 i + B 22 ){np+) - (B n - B 12 + B 21 - B 22 ){hpJ)\ 



2 

~ [(B n + Bi 2 - B 21 - B 22 ){np+) - (B n - B 12 - B 21 + B 22 ){hp. 



(B15) 



(B16) 



In Appendix C we prove 
B u + B 21 = l, B 12 + B 22 = l, 0<B^<1, K lk + K 2k = 0, K kl = K lk , 

which suggests the definitions 
B + = B n -B 22 , B_ = 1 - B n - B 22 , \B+\ < 1, |B_| < 1, |B + ±B_|<1. 

(B17) 

By taking into account these relations we obtain 

2 2 L 

KiBij(npj) - K iKu{n^i) = K + (hp + ) - K-[B+(hp + ) + B^(hp-)} 

i,j=l i=l 1=1 

L 

+ 5>-tfii("6) (B18) 
i=i 

so that the expression for the 1PI vertex function takes the form 

r°° L ( 1 d \ 
G(p\0(K:,n)) = J2 da]JP l (-—)e-^ m '-^f r (a;p;0 
r Jo l=l \2idZ,i) 

x expi |K + (np + ) - K_[B + (np + ) + B_(np_)] + ^ /^^(n^) j ||=o. 

(B19) 

The differentiation with respect to £ generates the scalar products PiPj, K-(np + ) as well as 
K_(np_) so that we are able to write 
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G(p\0(k; n)) = J2 dag r (a] PiPf, K_np)e 



in+{np+)— in- [B+ (hp+)+iB- (np_ )] 



(B20) 



where 



gr(a; PiPj ;K-np) = f[P t (^.-^ e^^^g^ 

xe i{K + (np + )-K_[B+(np + )+_B_ («P-)]+X)j=i |^ =Q 



(B21) 



According to the definition of the anomalous dimension ( |B3"|) , ([B4]) and (|B7|) , we perform 
the Fourier transform and find 



, d f d(np + ) f d(np 
7(«, k) = g— res 



■p ^ 



dg "J 2ix J 2tt 
x expi ^K+(np+) — K_[B + (np + ) + B_(np_)] — K r + (hp + ) — «/_(np_) j , 

(B22) 

where the i?'-operation acts onto the coefficient function g? only. For convenience, we 

(B23) 



introduce the variables 16 



K, — Kj 



ill 



and receive 



Jg-^resR' J Q dag v [or, PiPf. 



d d \ 
*a — ' 5 ( B + + W +) 5 ( B - + w - 

OW4- OVJ- I 



Because of the restrictions (|B17| ) of the variables B + and B- we conclude 
7(u> + , 7^ for \w±\ < 1, \w + ± W-\ < 1 



(B24) 



(B25) 
(B26) 



Note that the introduction of 7(u> + ,w7_) was suggested by the special variable dependence 
ofEq. flB24j ). 

Let us turn back now to QCD in axial gauge. As already mentioned, we have to apply a 
very complicated gluon propagator. It is in the spirit of the foregoing considerations to use 
the following a-representation 



00 POO 



D ab = Sab I da I d(3 



JO 



+ n v 







2 V 92if„ d2i^J 



(B27) 
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If we insert this propagator instead of the usual one into the a-representation of a diagram, 
then all investigations can be performed in a similar manner as before. However, it is more 
instructive to use instead of the (a, /^-representation ( |B27p the formal expression 



-5, 



ab 



da 



- ff- 



d 



xe ia(k 2 +ip)+i2£k 



+ n u 



d 



IT 



(B28) 



By writing now the a-representation of a graph, we have to include also the pre-factors of 
the gluon propagators 



d 



1 

- I TV 

2 1 m%i v 



i 



+ 



1 



^ 9 



'd2i£_ 



IT 



(B29) 



into the original product of the polynomials, so that instead of a polynomial there appears 
finally a rational function in the ^-derivative. This may introduce new singularities. How- 
ever, the expression in the exponential which is responsible for the support properties of the 
anomalous dimensions remains unchanged. This establishes the general support restrictions 
( [B25D for QCD too. 

Of course, for special operators there exist additional restrictions. The operator ( |B~T|) 
transforms under charge conjugation C in the following way (see Appendix D) 



C (O a (K+, k_; h)) C + = TO a (n + , -k_; h). 



(B30) 



If we apply the charge conjugation to the renormalization group equation (|B2[ ) then we get 
d 



H— O a {K + ,-K_-h) = J d 2 il (j(k,iI) -2 7 ^ (2) (k-kO) O a {n' + ,-K'_; 



n 



(B31) 



The substitution and the comparison with the definition of the operator 

( [BID and Eq. ([B2] ) leads to j(k+, k' + , k'_) = — k' + , —k'_), which is equivalent to 



7(U7+,U7_) = 7(-U7+,W_). 



(B32) 



APPENDIX C: PROPERTIES OF THE COEFFICIENTS 
OF THE a-REPRESENTATION 

In this Appendix we prove the relations ( |B16| ). We study the coefficients D(a), Bu(a) 
and Kki(a) of the a-representation of lPI-Feynman graphs with an operator insertion. We 
assume that the graph contains V vertices (the vertex corresponding to the operator insertion 
is the first vertex), N external and L internal lines. To each internal line there corresponds 
one a parameter. According to the momentum flow we define an orientation of the diagram. 
Instead of the coefficients Bn and appearing in the a-representation we investigate the 
coefficients 

ki (a) = B u (a) D (at) , k M (a) = K u (a) D (at) , (CI) 



35 



which are homogeneous polynomials with respect to all a and linear functions of each a [22 . 
They are given by topological formulas. D is the so-called chord set product sum 



Te{T} \1<£T J 

Here, {T} denotes a set of trees. A tree of the lPI-graph is generated by cutting (L — V) 
lines of the graph so that a connected graph remains which does not contain loops. The 
degree of homogeneity of D is L — V. We note further definitions. If we consider a special 
line I, then a tree containing this line is denoted by TJ. If we cut this line, then we divide the 
considered tree into two half-trees corresponding to two simple connected diagrams. The 
half-tree lying in the direction of the orientation of the line I is denoted by TJ + . 

After these explanations we write down the definition of the coefficients bu(a) with the 
degree of homogeneity (L — V) 



bu{a) = £ II a k 

Tie{T} \k<£T 



(C3) 



Here, the sum runs over all trees containing the line I whereas the considered half-tree has 
to contain the external vertex i. The set of these trees is denoted by {T/}| igT +. Since 
{^}|ieT+ c {^} an d a k > for each k, it follows from the definition of D(a) and bu(a) that 



< B u (et) 



biija) 
D(a) 



<1, 



(C4) 



which is the first relation (pi6| ). 
To prove the identity 

Bu(a) + B 2i (a) 



bu(a) + b 2 i(a) 
D(a) 



,V, 



(C5) 



we have to show that {T} is just the disjoint union of the sets {T\}^ ieT + and {T2\\ ieT +- Since 
T G {T} is connected, it has to contain the line l\ or l 2 (see Fig. |J). Therefore T is contained 
in the sets {7\}| ieT + or {T 2 }\ i€T +. Moreover {T} = {Ti}\ ieT + U {T 2 }^ ieT +. It remains to show 
that the sets {T±}\ ieT + and {T 2 }| igT + are disjoint. If the vertex % is contained in a tree T 
belonging to the sets {Ti}\ ieT + and {T 2 }\ ieT + then it is possible to find a path from the 
vertex i over the line l\ to the operator vertex and over the line l 2 back to the vertex i. But 
this is not possible because T cannot contain a loop. Therefore the set {T} is the disjoint 
union of the sets {Ti}\ ieT + and {T 2 }| igr + which proves Eq. (p5|). 
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FIG. 4. Topology of the lPI-Feynman graphs for the light-cone operators in light-cone gauge. 
The box symbolizes the connected (N + 2)-point function. The lines l\ and I2 are contained in a 
loop c. 



For the investigation of the properties of the Kki(a) we start from the definition 

fc(&0= E&Ms06= £ (n«*) (sw)) ■ ( C6 ) 

M=i T c e{T c } WT C / \jec J 

Here, T c denotes a pseudo tree (this means a tree to which we add a line in such a way that 

the new structure contains a loop c). The sum fe jec (±^)^ runs over all lines j which 
form the loop c. When the line j G c points along the orientation of the loop c then we 
have to write a plus sign, in the opposite case a minus sign. Obviously the differentiation 
with respect to £ provides the symmetry relation K^i{a) = Kik(a). Furthermore, the £1- 
and ^-dependent part of relation ( |C6| ) is determined by those T c in which l-y and I2 are 
contained in the loop c. Then, because of our conventions, the lines l\ and I2 have different 
orientations so that we can write 

fe(&i) = 6, 6, ■■■,&)■ (C7) 

The desired relation Kn(a) = —K 2 i(a) is then a simple consequence of Eq. ( |C7|) . 



APPENDIX D: CHARGE SYMMETRY 



We prove the following behaviour of the renormalized operators O a (K +1 k_) under charge 
(C)-conjugation: 

CO a (K + , n_)e iSl C + = TO a ( K+ , -«_). (Dl) 
Let us begin with the C-conjugation of the elementary fields. 

C^(x)C + = E c a ^(x), CA»(x)C + = -A»(x), 

P 

C^ i a (x)C+ = J2^)cp 1 a, (D2) 

$ 

where i and a, (3 denote the flavour and the spinor indices of the quark field if). The 

T 

transposed matrix of c := \c a p} is denoted by c. Using the relation 

c" 1 7 M c =\ (D3) 
we find for the charge conjugation of the operator : ip(x)'j^\ a U(x, y)ip(y) : 

C ^(x) lfl X a U(x } y)^(y):C + = E :C(x)(c _1 v) a ^A» [C{/(x,t/)C + g(2/): 

a/3 

= E -^)(^) 0a >^[CU(x i y)C + ]^: 

a/3 

= -^(y) 1 ,\ a U(y,x)^x):. (D4) 
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where we have used the anticommutation rule of the fermionic field as well as 

CU(x,y)C + = Pexp J" dz^CA^zjC^ = U(y,x). (D5) 

T 

The hermitean generators A a satisfy the relation A a = ±A a . Consequently, setting in Eq. 
x = K-\n, y = K 2 n, we find 

C :ip(Kih)'y f j i \ a U(Kin, K2n)if)(K 2 n) : C + = =F '■if){K 2 n) n f tl \ a U{K 2 n, K\n) if) (Kih) : 

(D6) 



or 



CO a (K + ,K^)C + = TO a (K + ,-K_), re±= (k 2 ±«i)/2. (D7) 
In QCD C-invariance leads to relation (|Dl|), 



APPENDIX E: LIMITS 

After some algebra we obtain the following results for the structures appearing in the 
BL-kernels: 

LIM- = l-z, (Ela) 



LIMF(x,y) = ^^, (Elb) 



UMF(x, y) In [ - ] = i±*L \ n ( z ) + i _ z (Elc) 

\yj l ~z 



LIM^ In f - ) = -i ((1 + z) ln(z) + l-z), (Eld) 



LIM (f(x, y) + J J In 2 ( ~J = I±^L ln 2 (^) + (1 - *) In(z), (Ele) 



LIM[F(x, y) - F(x, y)\ In ( - j In M - - J = 2^-^ ln(z) ln(l - z), (Elf) 



LIM[F(x, y) hx{y) hx{x) - F(x, y)U 2 (y) + F(x, y)Li 2 (x 
1 + z 2 /V_ 1 
~ 1 - z V 6 2 



' h : " % - ~ ln 2 (z) ] +(l + z) \n(z) + 2(1 - z), (Elg) 
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L1MF(x,y) 



x 



2Li 2 1 - - + ]n\y) - 2 ln(x) ln(y) + 2U 2 (y) - 2U 2 (x) + In 



2 i+^ u 

l + \z\\ 





z 




1 + 


z\ 



Li; 



+ 2Li 2 \l - - 

. 1 + N .+^ln a (|z|)-ln(|z|)ln(l + |^ 
+ 2(1 + |z|)1b(|*|)+4(1-|*I), 



1 



(Elh) 



LIM 



x 



ln(x)(l + ln(x) - 21n(x)) + 



\ x <- 




f] 




-> X 





0. 



(Eli) 



Whereas the determination of most of these expressions is very simple, those containing the 
Li2-functions are not so easy to handle. 

To determine the limit of G(x,y) in Eq. ( |5.8f ) we have to prove Eqs. ( |Elgj ) and ( |Elh| ). 
First, the expression 



lim r" 1 (2F[\n(y) ln(x) + U 2 (y) + Li 2 (x)] - 2F[ln(y) ln(x) + Li 2 {y) + Li 2 l 

T— >0 N 



X 



7T I X 

\n(y) ln(x) + Li 2 (y) + Li 2 (x) = — + Li 2 (x) - Li 2 (y) - ln(y) In = 

6 \v, 



ln(y) ln(x) + Li 2 (y) + Li 2 (x) 



7T 



X 



+ Li 2 (y) - Li 2 (x) - ln(x) In = 



(E2) 



may be simplified. Taking into account Li 2 (x) + Li 2 (l — x) — — ln(x)ln(l — x) + vr 2 /6, we 
get 



(E3) 



From this, together with [ 33 1 



1 



Li 2 ( r ) = r+-V + 



Li 2 ( 1 - - 



-Li 2 

I - iln 2 
z 2 



T — Z 



2 



Z — T 



T 



7T 



Z — T 



T 



71" 

if 



and the corresponding expansion of F and F, we obtain the result (|Elg ). 
Consider now the term within the second bracket of ( |5.81| ), 



x 



x 



9 - \9 1 - - sign(y) 



x 



(E4) 



2(F - F)Li 2 1 1 J +(F - F) ln 2 (y) - 2F\n(x) \n(y) 



+ 2FU 2 {y) - 2FU 2 (x) 
According to our procedure we have to calculate 



J X <r 


-> X 1 






\y<- 


J 



(E5) 
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E(z) = limr" 1 



X 



/ 'u . u) [ 2U, ( 1 - I + ln 2 (y) - 21n(x) Info) + 2Li 2 (y) - 2Li 2 (x) 



+ ln 2 (y) + 2Li 2 1-= - 



)H 


a; *r 




u 






-»• 37. 





P=f ,1/ 



T T 



Here, we apply Li 2 (x) + Li 2 (l/x) = — \ ln 2 (— x) — 7r 2 /6 to express 



Li 2 £ 



x 



X \ 7T 



In - In 1 - - + — 



-Li 2 



2/ 



l ln 2 



7T 
"6 



and the relation Li 2 (l — x) + Li 2 — ^ = — | ln 2 (x) to get 

1 



Li 2 1 - 



-Li 2 



x 



x - y 



In" 1 - 



x 



Then we obtain 



E(z) 



lim r 1 



2F ^Li 2 (y) -Li 2 (x) - Li 2 



y 



-2F Li 2 (x) -Li 2 (y) - Li 2 



x - y, 
f Li 2 



Li 2 



x 



x-y 



^ + ln(y — x) In ^ 



x 



+ In 



y-x 



yx 



\n[ y - 

x 



(E6) 



(E7) 



(E8) 



\X=-,V=-- 

(E9) 

From this transformed expression we are able to determine (|Elh|) using (|E4|) together with 

Li 2 (x + e) = Li 2 (x) - eln(l — x)/x H . (E10) 
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